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1. INTRODUQAO 

Os geradores, transform ado res e linhas de transmissao sao componentes 
essenciais em urn Sistema Eletrico de Potencia (SEP). Dessa forma, o 
monitoramento destes equipamentos deve ser constante para assegurar sua correta 
operagao, o fornecimento de energia e minimizar os custos decorrentes de uma 
eventual parada. Atualmente, alem dos tradicionais e diversos dispositivos de 
protegao, os equipamentos do Sistema Interligado Nacional (SIN) sao monitorados 
por Registradores Digitais de Perturbagao (RDP), conhecidos tambem como 
oscilografos digitais. Estes equipamentos monitoram constantemente grandezas 
predefinidas do SEP, como por exemplo, as tensoes, correntes e frequencia, bem 
como os estados de dispositivos de protegao, chaves seccionadoras e disjuntores. 
Caso ocorra alguma variagao nesses sinais que ultrapasse urn limiar pre- 
estabelecido, ou caso ocorra atuagao de fungoes de protegao ou chaveamentos, o 
equipamento cria urn registro que pode ser acessado posteriormente (VARELA et 
al., 2009). 

Urn registro de oscilografia e composto por urn conjunto de dados digitais 
dispostos em urn formato padronizado e nomeado de COMTRADE pelo Institute of 
Electrical and Electronic Engineers (IEEE) (IEEE STD, 2013). O padrao IEEE 
COMTRADE (Common Format for Transient Data Exchange) define urn formato 
comum para os arquivos de dados digitais e assim urn meio para intercambio de 
informagoes. Cada registro COMTRADE possui urn conjunto de ate quatro arquivos 
associados. Cada urn dos quatro arquivos contem diferentes informagoes (SANTOS, 
2011). 

Em todo o processo de analise dos registros de oscilografias, o tempo de 
coleta dos dados e filtragem dos registros relevantes pode durar algumas horas 
(GIOVANINI et al., 2008). Somado a isto esta o fato de que diversas ocorrencias 
podem ser registradas por dia em uma empresa de geragao ou transmissao de 
energia eletrica, sendo a maioria correspondente a procedimentos operativos 
normais, como energizagao e desenergizagao de geradores. 

Os RDPs sao ajustados com tal sensibilidade que garanta o registro de 
qualquer anormalidade no sistema, mesmo que esta nao seja suficiente para 
produzir urn blecaute de grandes proporgoes. Assim a quantidade de informagoes 
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registradas pelos RDPs e muito grande, produzindo uma demanda de trabalho para 
analise manual que por vezes a torna inviavel (MORETO e ROLIM, 2010). 

No atual modelo de operagao do SIN o ONS desempenha a fungao de 
ditar as regras as quais os agentes do setor devem se submeter. Dentre estas 
regras, que sao denominadas Procedimentos de Rede, encontra-se o Submodulo 
1 1 .6 que trata dos registros de perturbagoes. Este submodulo atribui aos agentes a 
responsabilidade de implantar urn sistema de coleta de dados provenientes dos 
RDPs (ONS, 2009). Como os RDPs sao dispositivos digitals essa coleta se da 
atraves de uma rede de comunicagao que permite acesso remoto aos dados e 
interliga os RDPs a urn servidor remoto que centraliza as informagoes em urn unico 
local, como exemplifica a Figura 1 . 



Servidor 




_U_sina_l 

Figura 1 : Rede de coleta de registros. 
Fonte: (MORETO e ROLIM, 2010). 

E perceptivel que essa rede de RDPs concentra e disponibiliza uma 
grande quantidade de dados aos analistas, gerando uma alta demanda de trabalho. 
Considerando que os oscilografos ou RDPs produzem registros dos sinais 
monitorados em formato digital torna-se possfvel o emprego de metodos 
computacionais para promover a analise destes registros. O proprio submodulo 1 1 .6 
do ONS recomenda aos agentes: 



"E permitida a utilizacao de filtros automaticos, instalados no agente de 
transmissao, para descarte de registros ocasionados por pallidas 
programadas dos RDP, para a verificacao da condigao de operagao e 
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calibragao desses registradores ou para a eliminagao de partidas espurias, 
tais como as decorrentes de manobras de secionadoras." (ONS, 2009). 

A rapida analise dos eventos, pela pos-operagao, e de suma importancia 
para diagnosticar a sua causa, permitindo agir sobre um possivel defeito, 
minimizando a indisponibilidade de equipamentos ao SEP. Como o tempo 
despendido nas analises e crucial para a diminuigao dos custos decorrentes da 
parada de um equipamento, e evidente a necessidade de uma metodologia capaz 
de gerar um diagnostico da perturbagao de forma automatica, apontando para uma 
possivel causa do problema e possibilitando priorizar a atengao do engenheiro 
especialista aos casos mais importantes. Dessa forma, o processo de analise de 
ocorrencias se torna mais agil e reduz-se o tempo de indisponibilidade de 
equipamentos. (MORETO e ROLIM, 2010) (DAVILA, 2010) (PEREIRA e SPRITZER, 
2007). 

O desenvolvimento de um sistema de analise automatica de oscilografias 
deve levar em consideragao o objetivo para o qual a analise sera direcionada. A 
finalidade da analise ira determinar os dados (insumos) do registro de oscilografia a 
serem utilizados (registros de curta duragao, registros fasorial, sequencias de 
eventos). Antes de um sistema conseguir realizar uma analise faz-se necessario 
ainda uma etapa de extragao de dados dos registros. Esta etapa utiliza ferramentas 
matematicas para transformar as variagoes numericas nos dados registrados em 
informagoes que possam ser comparadas com padroes de causas ja conhecidas, 
conseguindo assim determinar o que gerou a perturbagao. 

Nesse contexto alguns metodos tern sido propostos, utilizando diferentes 
tecnicas e ferramentas computacionais, para realizar a analise de forma automatica 
dos registros dos RDPs. Lenz (2009) descreve o desenvolvimento de uma 
ferramenta computacional para avaliar o desempenho de reles de protegao usando 
os dados dos RDPs e das Sequencias de Eventos (SOE), produzidas pelos sistemas 
de supervisao e controle. Costa (2010) propos uma tecnica de diagnostico de 
disturbios transitorios no sistema de transmissao baseada na transformada wavelet 
para aplicagao integrada ao desenvolvimento de novos RDPs. Moreto (2011) 
desenvolveu um sistema de analise automatica para registros digitais de 
perturbagoes em unidades geradoras usando limiares exatos para classificar as 
causas dos eventos. 
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O foco deste Estudo Especial e assimilar conceitos matematicos da 
transformada wavelet, bem como levantar o estado da arte quanto a aplicagao deste 
metodo no processamento digital de sinais, estudar as ferramentas computacionais 
para implementagao do metodo matematico "transformada wavelet" manipulando 
dados dos registros de oscilografias e extrair as caractensticas necessarias dos 
sinais; 

O capitulo 2 deste estudo especial faz uma revisao sobre oscilografia 
explicitando o formato dos dados registrados e as possiveis aplicagoes para a 
analise destes dados. Na sequencia, o capitulo 3, aborda a transformada wavelet, 
ferramenta matematica amplamente utilizada na analise de oscilografias. 



11 



2. OSCILOGRAFIA 

O termo "oscilografia" refere-se aos registros produzidos pelos 
oscilografos. Estes sao equipamentos desenvolvidos para registrar oscilagoes 
eletricas e o nome vem de "oscilar" + "grafo", ou seja, escrita de oscilagoes 
(WEISZFLOG, 2004). Originalmente estes equipamentos eram dotados de agulhas 
registradoras conectadas a bobinas eletromagneticas que transformavam oscilagoes 
de corrente eletrica em oscilagoes mecanicas na ponta de uma agulha, que por sua 
vez, literalmente, escreviam o registro da oscilagao. 

Com a evolugao tecnologica os antigos oscilografos se modernizaram e 
passaram a ser chamados de Registradores Digitais de Perturbagoes, ou pela sigla 
RDPs, que atualmente se baseiam em tecnologia digital e mantem os registros em 
memoria semicondutora. Porem, os profissionais do meio tecnico e academico se 
mantem fieis ao uso do termo oscilografia para se referir aos arquivos digitais 
produzidos pelos RDPs. A Figura 2 mostra urn RDP moderno da Reason Tecnologia 
S.A. 




Figura 2: Registrador Digital de Perturbagoes Multifuncoes RPV-310. 
Fonte: (REASON, 2010). 
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Este modelo tem sido largamente empregado pelas concessionarias de 
energia e possui as seguintes especificagoes tecnicas (REASON, 2010): 

• Ate 64 entradas analogicas (tensao, corrente e transdutores DC); 

• Ate 384 entradas digitais; 

• Ate 256 grandezas binarias obtidas de mensagens GOOSE segundo IEC 
61850; 

• Conversores analogico-digitais de 1 6 bits; 

• Taxa de amostragem de 256 pontos-por-ciclo; 

• Gravagao de registros de forma de onda com 256, 1 28 ou 64 pontos-por-ciclo; 

• Registros de longa duragao continua e por trigger; 

• Aquisigao sincronizada a referenda externa de tempo padrao IRIG-B; 

• Medigao e registro de harmonicas ate 50 a ordem; 

• Sinalizagao de falha do equipamento. 

No formato digital existem diversas maneiras de gravar informagoes e 
diferentes fabricantes de equipamentos tendem a utilizar diferentes tecnicas. Porem 
isto dificulta o acesso e manipulagao destas informagoes pelos usuarios. Para 
facilitar a troca de informagoes dos registros e compatibilidade entre fabricantes o 
ONS adota o padrao IEEE COMTRADE para padronizar os registros dos RDPs 
instalados no sistema eletrico brasileiro (ONS, 2009). 

O padrao COMTRADE, determina que as informagoes dos RDPs sejam 
divididas e gravadas em quatro arquivos de dados. Os quatro arquivos devem 
possuir o mesmo nome apenas diferindo na extensao, a qual indica o tipo de 
conteudo do arquivo. Os arquivos com suas extensoes sao os seguintes (IEEE 
STD., 1999): 

• Arquivo de cabegalho (xxxxxxxx.HDR) - e urn arquivo de texto no formato 
ASCII. As informagoes contidas neste arquivo podem ser impressas e 
manipuladas diretamente pelo usuario que pode entender melhor as 
condigoes apresentadas. Nao ha urn padrao a ser seguido, o IEEE apenas 
sugere alguns tipos de informagoes possiveis para serem descritas em urn 
arquivo de cabegalho, as quais irao auxiliar na interpretagao dos dados: 
descrigao do sistema de energia, nome da estagao, identificagao da linha, 
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transformador, reator, capacitor, ou disjuntor que e monitorado por cada 
canal, tensao nominal de enrolamento de transformadores, especialmente a 
potencia e a corrente, a sequencia de fase de entrada dentre outras 
informagoes (SANTOS, 2011); 

• Arquivo de configuragao (xxxxxxxx.CFG) - e um arquivo de texto no formato 
ASCII que e lido por softwares especializados, por isso, e escrito em formato 
especifico. O arquivo e dividido em linhas e as linhas em campos. Virgulas 
sao utilizadas como delimitador de campo, mesmo que o respectivo dado nao 
seja especificado permitindo que o comprimento de cada campo seja variavel. 
A informagao em cada linha do arquivo deve ser listada na ordem exata 
requerida. Qualquer desvio do padrao estipulado anula todo o conjunto de 
arquivos. Contem as informagoes necessarias para interpretagao do arquivo 
de dados (.DAT) e por isso nao e opcional. As seguintes informagoes devem 
estar contidas no arquivo de configuragao (IEEE STD., 1999): 

o Nome da estagao, identificagao do dispositivo de monitoragao, ano de 
revisao do padrao COMTRADE; 

o Numero e tipos dos canais; 

o Nome dos canais, unidades, a fatores de conversao; 

o Frequencia da linha; 

o Taxa de amostragem; 

o Data e tempo do primeiro dado; 

o Data e tempo do ponto de acionamento da falta (trigger); 

o Tipo do arquivo de dados; 

o Fator de multiplicagao da estampa de tempo. 
A Figura 3 mostra um exemplo de arquivo de configuragao, onde pode ser 
visto na primeira linha a identificagao da subestagao, do dispositivo de 
monitoragao e o ano do padrao COMTRADE. A segunda linha especifica nove 
canais de entrada sendo todos analogicos, e nas linhas seguintes estao os 
nomes dos canais com suas unidades e fatores de conversao. Na sequencia e 
especificada a frequencia, a taxa de amostragem, data do primeiro dado e do 
ponto de trigger, tipo do arquivo e fator de multiplicagao da estampa de tempo. 
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SE,RDP, 1999 
9, 9A, OD 

1, VA, , ,V, 17 .44, 0, 0, -327 68, 327 67, 1, 1, S 

2, VB, , ,V, 17 .44, 0, 0, -327 68, 327 67, 1, 1, S 

3, VC,,,V, 17. 44, 0,0, -32768, 32767, 1,1, S 

4, Reator_IA, ,, A, 0.221, 0,0, -32768, 32767, 1,1, S 

5, Reator_IB, ,, A, 0.221, 0,0, -32768, 32767, 1,1, S 

6, Reator_IN, ,, A, 0. 05524, 0,0, -32768, 32767, 1,1, S 

7, LT la, A, , A, 0.023187, -1.28851, 0,-32767, 32767, 1200, 5, S 

8, LT lb, B,, A, 0. 0637747, -145. 13, 0,-32767, 32767, 1200, 5, S 

9, LT In, N, , A, 0 . 08 64 938, 227 . 923 , 0 , -327 67 , 327 67 , 12 0 0 , 5 , S 
60 

1 

1440,1968 

21/02/2007,22:35:56.8 95000 
21/02/2007, 22: 35: 57. 14667 4 
ASCII 

1 

Figura 3: Exemplos de registro de oscilografia: arquivo de configuragao (xfg) 
Fonte: Autoria propria. 



• Arquivo de dados (xxxxxxxx.DAT) - contem o numero da amostra, a estampa 
de tempo e o valor de cada amostra dos canais de entrada do registro 
gravado. Os dados sao dispostos conforme descreve o arquivo de 
configuragao que tambem especifica os fatores de conversao para converter 
os valores apresentados em dados concretos, pois os numeros armazenados 
representam apenas uma escala dos valores reais. Alem dos dados das 
entradas analogicas, pode-se registrar entradas digitais. Nese caso, o estado 
da entrada e associado a mveis logicos (SANTOS, 201 1). 

O arquivo de dados pode ser escrito em formato ASCII ou em formato 
binario e isto e especificado em urn dos campos do arquivo de configuragao. 
Para o formato ASCII o delimitador de campos e a virgula enquanto em 
arquivos binarios usa-se o tamanho de cada informagao como delimitador. 

No formato ASCII, o arquivo e dividido em linhas e colunas. O numero de 
linhas varia de acordo com o tamanho da amostra e afeta diretamente o 
tamanho total do arquivo. O numero de colunas depende da quantidade de 
canais analisados e tambem afeta o tamanho total do arquivo. A seguir e 
mostrada a sequencia de informagoes necessarias (IEEE STD., 1999): 

o A primeira coluna contem o numero da amostra; 

o A segunda coluna e a estampa de tempo; 

o A terceira coluna representa o conjunto de dados dos canais 
analogicos; 

o As proximas colunas representa o conjunto de dados dos canais 
digitais. 
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Nao pode haver quebra de linha entre os dados referentes ao mesmo 
tempo de amostragem. O arquivo deve terminar com urn caractere ASCII de 
EOF, End Of File, referente ao termino de arquivo. 

A Figura 4 mostra o conteudo do arquivo .DAT para o mesmo registro da 
Figura 3. Verifica-se que o arquivo e composto por onze colunas, separadas 
por virgulas, conforme o formato ASCII especificado no arquivo de 
configuragao. A primeira coluna especifica o numero da amostra, a segunda a 
estampa de tempo e as nove seguintes o valor da amostra para o canal 
correspondente. Cada linha corresponde a uma nova amostra para os 
mesmos canais. 



0000000001, 0000000000, 00 6343, 004197, -11227, 000216, -00234, 000026, 000055, 002275, -02 635 
0000000002, 00000006 94, 004218, 00 6300, -11082, 000259, -00191, 000025, 000055, 002275, -02 635 
0000000003, 0000001388, 001707, 008611, -10502, 000252, -00155, 000016, 000055, 002275, -02 635 
0000000004, 0000002083, -01557, 010612, -08871, 000264, -00084, 000024, 000055, 002275, -02 635 

0000000005, 0000002777, -04251, 011157, -0 6455, 000247, -00027, 000019, 000055, 002275, -02 635 

0000000006, 0000003472, -0 6356, 011005, -03 97 8, 000206, 000051, 000023, 000055, 002275, -02 635 
0000000007, 0000004166, -08615, 010421, -01432, 000151, 000110, 000019, 000055, 002275, -02 635 
0000000008, 0000004861, -10590, 008 652, 001612, 0000 92, 0001 67, 000024, 000055, 002275, -02 635 
000000000 9, 0000005555,-11324, 00634 9, 004403, 000035, 000222, 000024, 000055, 002275,-02635 
0000000010, 000000 6250, -1133 9, 004231, 006469, -00039, 000255, 000027, 000055, 002275, -02 635 

0000000011, 0000006944, -10 675, 001713, 00857 9, -00098, 000268, 000019, 000055, 002275, -02 635 

0000000012, 0000007 638, -08772, -014 62, 010547, -00166, 000257, 000029, 000055, 002275, -02 635 

0000000013, 0000008333, -0 6353, -04197, 011230, -00210, 000239, 000022, 000055, 002275, -02 635 
0000000014, 0000009027, -04229, -0 62 94, 011084, -00242, 000201, 000019, 000055, 002275, -02 635 
0000000015, 000000 9722, -01707, -08608, 010487, -00260, 000152, 000027, 000055, 002275, -02635 
0000000016, 0000010416, 001551, -10610, 008867, -00253, 000101, 00002 9, 000055, 002275, -02635 
0000000017, 0000011111, 004244, -11159, 006467, -00230, 000025, 000025, 000055, 002275, -02635 
0000000018,0000011805,006348,-11007,003985,-00174,-00030,000024,000055,002275,-02635 

Figura 4: Exemplos de registro de oscilografia: arquivo de dados (.dat) 
Fonte: Autoria propria. 



• Arquivo de informagoes (xxxxxxxx.INF) - e urn arquivo opcional em formato 
ASCII que contem informagoes extras que a fonte de dados deseja tornar 
dispomvel para os usuarios. Este arquivo tambem pode conter informagoes 
privadas, so interpretaveis pelo fabricante (IEEE STD., 1999). 

Alem das informagoes contidas nos arquivos COMTRADE, muitos 
fabricantes de equipamentos implementaram uma funcionalidade a mais, o registro 
de longa duragao ou simplesmente registro fasorial. Esse registro apresenta uma 
duragao da ordem de minutos e registra apenas informagoes de modulo e angulo 
das tensoes e correntes amostradas em uma taxa menor. Essas informagoes sao 
uteis para analises de fluxo de potencia e oscilagoes de baixa frequencia, como as 
de origem mecanica (MORETO, 201 1 ). 
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Como esse tipo de registro nao estava definido no padrao COMTRADE 
definido na norma IEEE C37.1 1 1 -1999, cada fabricante de RDP criou seu proprio 
padrao, acrescentando arquivos ao registro. A nova revisao do padrao IEEE 
C37.1 11-2013 define um esquema para acrescentar as informagoes desse registro 
no arquivo da dados (.DAT) adicionando dois canais analogicos, um para o modulo e 
outro para a fase, de cada sinal que se deseja char um registro fasorial (IEEE STD, 
2013). 

O conteudo dos registros no padrao COMTRADE pode assim ser lido por 
softwares especializados, que utilizam as informagoes do arquivo de configuragao, 
para interpretar os arquivos de dados e construir graficos. Assim a informagao torna- 
se amigavel aos olhos de um analista, que consegue interpretar as oscilagoes 
registradas e determinar uma causa para a mesma. 

Ao transformar os dados dos registros de oscilografia para informagoes, 
podemos classifica-las em tres tipos principals (MORETO e ROLIM, 2010): 

• Registro de curta duragao - consiste no grafico produzido atraves dos dados 
dos canais analogicos de tensao e corrente monitorados, com uma taxa de 
aquisigao alta, da ordem de 256 pontos-por-ciclo (REASON, 2010). Em 
fungao da alta taxa de aquisigao sao armazenados apenas alguns segundos. 
Este registro tambem e chamado de registro de forma de onda, pois permite 
visualizar a forma de onda do canal monitorada. A Figura 5 mostra um 
exemplo deste tipo de informagao para tres canais de tensao monitorados, 
sendo um canal para cada fase. 



V-GER1 [kv] 

15 r r r r r r 




-15 t t t t t t t t t t 1 

2 2.01 2.02 2.03 2.04 2.05 2.06 2.07 2.08 2.09 2.1 

Tempo t[s] 

Figura 5: Exemplos de registro de curta duragao. 

Fonte: Autoria propria. 
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Registro de longa duragao - consiste no registro dos fasores de um canal, ou 
seja, apenas o modulo e o angulo do sinal monitorado sao registrados, com 
uma taxa de aquisigao bem menor que o registro de forma de onda, 
permitindo-se assim registrar ate alguns minutos de informagoes. A Figura 6 
mostra um exemplo de registro fasorial, como tambem e conhecido. 
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Figura 6: Exemplos de registro de longa duragao. 
Fonte: Autoria propria. 



Sequencia de eventos - O monitoramento dos canais digitais permite criar 
uma sequencia cronologica dos eventos ocorridos, como atuagao de 
protegao, abertura de disjuntor, atuagao de protegoes secundarias ou de 
retaguarda entre outros. A Figura 7 mostra um exemplo de sequencia de 
eventos em um registro de oscilografia. 
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: jVA04_87L86BFJ4( : 



-0.3 



Figura 7: Exemplos de registro de sequencia de eventos. 
Fonte: Autoria propria. 
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As leituras dos registros a que se referem as Figura 5 e Figura 6 foram 
realizadas em rotina desenvolvida no MATLAB® 1 e para a Figura 7 foi realizada no 
software SIGRA 4 - Fault Record Evaluation® 2 . 

2.2 APLICAQOES PARA A OSCILOGRAFIA 

A analise das ocorrencias e perturbagoes no SIN implica diretamente na 
seguranga operacional do sistema e esta relacionada com a qualidade dos servigos 
de energia eletrica. Por qualidade dos servigos entenda-se continuidade, 
confiabilidade, tempo restabelecimento em casos de blecaute e outros. Essa 
responsabilidade e dos agentes de operagao do setor eletrico e do ONS, de acordo 
com o definido pela a Resolugao Normativa Aneel n Q 115, de 29 de novembro de 
2004 (COSTA, 2010). 

As analises de ocorrencias e perturbagoes proveem subsidios 
fundamentals para identificar as causas e consequencias, para avaliar e especificar 
equipamentos de protegao, como reles e disjuntores, para o diagnostico de 
equipamentos e classificagao dos disturbios, entre outras aplicagoes. Desta forma, 
medidas corretivas e preventivas podem ser tomadas, resultando em melhorias no 
desempenho do sistema eletrico e da qualidade dos servigos publicos de energia 
eletrica (PEREZ e RINCON, 2012). 

As oscilografias constituem o principal insumo para as analises de 
ocorrencias e perturbagoes, sendo o seu registro considerado pelo ONS como 
documento oficial para o Relatorio de Analise de Perturbagao (RAP) (ONS, 2009). 
Atraves da analise de oscilografias e possivel extrair informagoes uteis para diversas 
aplicagoes. Os topicos seguintes fazem uma breve descrigao sobre alguns campos 
de pesquisa aonde as oscilografias verm sendo empregadas. 



1 MATLAB®: ambiente de desenvolvimento de algoritmos e analise de dados da MathWorks Inc. 

2 SIGRA 4 - Fault Record Evaluetion®: software de visualizacao de oscilografias da Siemens AG. 
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2.2.1 Analise da qualidade de energia 

A Agenda Nacional de Energia Eletrica (Aneel) define atraves dos 
Procedimentos de Distribuigao de Energia Eletrica no Sistema Eletrico Nacional 
(PRODIST) os criterios de qualidade rmnimos que os consumidores de energia 
eletrica tern direito. Especificamente no modulo 8 do PRODIST estao caracterizados 
os fenomenos de Qualidade da Energia Eletrica (QEE), os criterios de amostragem, 
os valores de referenda e os procedimentos relativos a qualidade do produto 
(ANEEL, 2008). 

Conforme o PRODIST sao considerados aspectos da QEE em regime 
permanente ou transitorio: 

• Tensao em regime permanente; 

• Fator de potencia; 

• Harmonicos; 

• Desequilibrio de tensao; 

• Flutuagao de tensao; 

• Variagoes de tensao de curta duragao; 

• Variagao de frequencia. 

As informagoes de todos os fenomenos acima citados podem ser 
extraidas de uma oscilografia de curta duragao. Fica assim explicito a aplicagao das 
oscilografias na analise da QEE, seja por parte dos consumidores ou dos 
fornecedores de energia. Kezunovic (2001) propos a viabilidade da automatizagao 
da analise de disturbios de qualidade de energia com uso de tecnicas de 
processamento de sinais aliadas a ferramentas de inteligencia artificial. A extragao 
de propriedades estatisticas dos sinais com uso de transformada de Fourier 
permitiria avaliar os parametros de qualidade de energia (KEZUNOVIC e LIAO, 
2001). Costa (2011) apresentou uma ferramenta para baseada em wavelet para 
detecgao de disturbios da qualidade de energia. As caractensticas da perturbagao 
podem ser identificadas por meio dos seus coeficientes de wavelet (COSTA et al., 
2011). 



20 



2.2.2 Diagnostico de equipamentos 

Monitorar o funcionamento de dispositivos de manobra, como disjuntores, 
ou de equipamentos como transformadores de instrumentos (transformador de 
corrente, transformador de potencial, para-raios) e uma constante necessidade no 
departamento de manutengao dos agentes do setor eletrico. 

Atraves da analise de uma oscilografia e possivel, por exemplo, 
determinar se o tempo de abertura de urn disjuntor encontra-se dentro do 
especificado pelo fabricante. Este tipo de equipamento mecanico de manobra sofre 
desgaste natural de seus componentes e definir o momento certo para realizar a 
manutengao preventiva e crucial para o sistema eletrico. Tambem e possivel 
determinar tempos de restabelecimento de urn para-raios apos a ocorrencia de uma 
descarga atmosferica e verificar se o mesmo nao apresenta corrente de fuga acima 
da nominal. Diversas outras caractensticas de equipamentos podem ser 
identificadas em uma oscilografia, como nivel de saturagao em transformadores de 
corrente, corrente de inrush de transformadores, reatores e banco de capacitores 
(MORETOe ROLIM, 2010). 

Urn sistema de monitoramento de disjuntor, para auxiliar na manutengao 
preditiva, foi proposto por Silva (2005). A transformada de Fourier foi usada para 
acompanhar a componente de frequencia fundamental das correntes e tensoes e a 
transformada wavelet para detectar o instante de interrupgao da corrente e assim 
determinar o tempo de abertura do disjuntor (SILVA, 2005). 

Lachman (2007) propos e implementou uma rotina computacional 
aplicando filtro de Fourier sobre registros de oscilografia para monitorar o 
funcionamento de Transformadores de Potencial Capacitivos (TPCs) instalados em 
subestagoes em operagao (LACHMAN, 2007). 



21 



2.2.3 Diagnostico dos sistemas de protegao 

Os requisitos basicos para um sistema de protegao sao (ANDERSON, 

1998): 

• Confiabilidade - habilidade da protegao em atuar corretamente quando 
necessario e evitar atuagoes desnecessarias; 

• Seletividade - capacidade da protegao em selecionar e interromper apenas os 
disjuntores mais proximos da falta; 

• Velocidade - rapidez com que a protegao localiza e elimina um defeito; 

• Sensibilidade - capacidade da protegao em atuar sob rmnimas condigoes de 
falta dentro da sua zona de protegao, sem perder a confiabilidade. 

Essas caractensticas dos sistemas de protegao podem ser monitoradas 
atraves das oscilografias. A velocidade pode ser verificada pelos tempos de atuagao. 
A seletividade pela eliminagao do defeito interrompendo somente os circuitos 
envolvidos na falta. A sensibilidade pela verificagao dos limiares de atuagao pre- 
ajustados nos reles de protegao e a confiabilidade pela verificagao da real 
necessidade da atuagao quando esta ocorrer (CARDOSO, 2013). 

Lens (2009) desenvolveu um programa computacional que relaciona os 
dados extraidos das oscilografias com as sequencias de eventos dos sistemas de 
supervisao. A ferramenta proposta avalia o desempenho dos reles de protegao apos 
uma ocorrencia apoiando o especialista em analise de protegao (LENZ, 2009). 

2.2.4 Localizagao de faltas 

Os sistemas eletricos devem oferecer alto grau de continuidade no 
fornecimento de energia. Porem, interrupgoes nos servigos de energia eletrica 
podem ser provocadas pela ocorrencia de diferentes tipos de fenomenos 
eletromagneticos que podem afetar distintas classes de consumidores. Perturbagoes 
podem ocorrer em diversos componentes de um SEP, dentre os quais se destacam 
as Linhas de Transmissao (LTs) como o elemento mais susceptivel, especialmente 
devido as suas dimensoes fisicas, complexidade funcional e ao ambiente em que se 
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encontram, apresentando uma maior dificuldade para manutengao e monitoramento 
(SILVA, OLESKOVICZ e COURY, 2005). 

As decisoes mais importantes sobre a operagao do SEP sao tomadas nos 
Centros de Operagao dos Sistemas (COS), onde em tempo real informagoes sobre o 
estado dos equipamentos sao recebidas por operadores humanos. Estes sao o elo 
final de uma prolongada cadeia de operagao e tern como fungao realizar rapidas 
analises de ocorrencias, buscando minimizar o tempo de restabelecimento das 
condigoes normals dos servigos de energia eletrica. Neste contexto ferramentas que 
indiquem a localizagao, mesmo que aproximada, de uma falta em uma LT, aceleram 
os procedimentos de manutengao. Equipes de campo especializadas podem ser 
acionadas e encaminhadas diretamente ao local da perturbagao no sistema, 
reduzindo o tempo de deslocamento (FARIA et al., 2009). 

Os dados registrados pelos RDPs, aliados ao conhecimento de 
caractensticas de uma LT, permitem a implementagao de tecnicas para calculo e 
determinagao da provavel localizagao da falta. O trabalho de Silva (2005) 
apresentou uma aplicagao da transformada wavelet para analisar os transitorios de 
alta frequencia em urn sistema de transmissao gerados por uma situagao de falta, 
tendo como principal objetivo, a determinagao do ponto de ocorrencia da mesma 
(SILVA, OLESKOVICZ e COURY, 2005). 

Os mais modernos RDPs ja tern disponiveis para o usuario fungoes de 
localizagao de falta baseadas na teoria de ondas viajantes. Em fungao do tempo de 
propagagao do sinal entre o ponto de falta e os terminals da linha, a localizagao da 
falta e estimada. Este e o caso do RPV-310 da Reason Tecnologia S.A. conforme 
descreve seu manual no topico sobre suas caractensticas: 

"Preparado para receber o modulo de localizagao de faltas por ondas 
viajantes. Tecnologia que permite a localizagao de faltas com a precisao de 
300 metros em linhas de transmissao / distribuicao de energia, 
independente do comprimento da linha" (REASON, 201 0). 

2.2.5 Detecgao e classificagao de faltas 

Dentre todas as possiveis aplicagoes para as oscilografias, uma delas se 
destaca. E a busca por determinar a causa ou tipo de falta que ocorreu no SEP, ou 
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simplesmente a classificagao da falta. A principal motivagao para esta aplicagao se 
deve as caractensticas dos SEP, que estao sujeitos a uma grande gama de defeitos 
de dificil determinagao pelas equipes de operagao em tempo real dos COS. Assim 
quando perturbagoes ocorrem e causam desligamento de fungoes do sistema, como 
geradores ou LTs, e suas causas nao sao facilmente detectadas, cabe aos analistas 
de oscilografia indicar ao menos uma provavel origem do defeito, para que as 
equipes de campo tenham urn ponto de partida no processo de manutengao. 

Esse processo pos-operagao de analise de oscilografias, geralmente 
dispoem de uma rede de coleta de oscilografias, que pode acessar todos os RDPs 
de urn agente do SIN, conforme mostrado na Figura 1 (MORETO e ROLIM, 2010). 
Esse processo centralizado de analise de oscilografias pode ser dividido em duas 
etapas, ambas com sua importancia, sendo elas explicitadas nos dois topicos 
seguintes. 

2.2.5.1 Detecgao de faltas 

A logica de funcionamento dos RDPs consiste no ajuste de condigoes de 
trigger. Estas condigoes quando atendidas dao initio ao processo de gravagao de 
urn registro de oscilografia. Os triggers nada mais sao do que limiares fixos para os 
canais analogicos monitorados e/ou mudanga de estado para os canais digitais. Por 
exemplo, pode-se configurar que se o valor de urn canal de tensao chegar a 0,9 pu 
dispare o processo de gravagao. Em urn canal digital o trigger pode representar a 
atuagao de uma protegao. Essa logica e muito eficiente e permite uma otima 
correlagao entre os dados de diferentes RDPs, por exemplo, instalados nos dois 
extremos de uma LT. (MORETO e ROLIM, 2010). 

O inconveniente e que uma unica perturbagao pode gerar o disparo de 
diversos RDPs no SEP, gerando registros em todos eles. Muitas ocorrencias que 
estao relacionadas apenas com fenomenos de QEE podem atingir os limiares de 
ajuste de trigger e produzir registros de oscilografia nao relacionados a uma falta. 
Tambem ocorre disparo de trigger devidos a manobras no SEP ou ensaios em 
equipamentos durante manutengoes preventivas. 

A etapa de analise de oscilografias com objetivo de detecgao de faltas 
consiste em uma analise menos aprofundada dos registros, que tern por objetivo 
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principal fazer uma triagem dos registros. Apos a triagem dois grupos de registros 
sao separados, um dos que nao representam faltas reais no sistema e outro que 
contem informagoes de faltas no sistema. Assim os registros relacionados a 
perturbagoes no SEP serao objeto de estudo mais minucioso enquanto os demais 
podem ser arquivados. O trabalho de Moreto (2011) buscou desenvolver uma 
metodologia capaz de arquivar automaticamente os registros de ocorrencias que 
nao necessitam de analise detalhada (MORETO, 2011). A Companhia Hidro Eletrica 
do Sao Francisco (Chesf) mostrou sua contribuigao com o desenvolvimento de um 
software responsavel pela triagem dos registros oscilograficos de seu sistema de 
transmissao de energia (BRITO, SOUZA e COSTA, 2010). 

2.2.5.2 Classificagao das faltas 

Na etapa posterior a triagem dos registros, busca-se fazer a classificagao 
das faltas. Esta analise tipicamente faz uso dos registros de curta-duragao, ou 
registros de forma de onda, que possuem alta taxa de amostragem. Faz-se 
necessario recorrer a tecnicas de processamento de sinais mais avangadas para 
auxiliar nessa tarefa, com o objetivo de automatizar o diagnostico sobre os tipos de 
disturbio presentes nos sinais registrados. A transformada wavelet constitui-se em 
uma ferramenta matematica bastante eficaz na analise de sinais de tensao ou 
corrente, obtidos durante a ocorrencia de disturbios na rede eletrica (CANDIDO, 
2008). Tambem se pode estimar o tempo/instantes das mudangas nos parametros 
do sinal registrado, segmentando-o em pre-falta, durante a condigao de falta e pos- 
falta. Apos a segmentagao do sinal em secgoes, o processamento e analise do sinal 
podem ser realizados nestes segmentos (UKIL e ZIVANOVIC, 2007). 

A classificagao pode indicar o tipo do defeito a que o SEP foi submetido. 
Para Meneses e Rocha (2009) as faltas em LTs podem ser classificadas como: 
curto-circuito fase-fase, trifasico ou fase-terra, rompimento de condutor, queimada, 
descarga atmosferica e outros (MENESES e ROCHA, 2009). O trabalho de Costa 
(2010) propos um metodo de diagnostico de disturbios transitorios baseado na 
analise em tempo real da energia dos coeficientes wavelets das tensoes e correntes 
em sistemas eletricos de transmissao de energia. O metodo e capaz de detectar os 
principals disturbios com transitorios, tais como faltas, faltas evolutivas, 
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afundamentos de tensao e transitorios devido as manobras de chaveamento em 
sistemas (COSTA, 2010). Ja o trabalho de Moreto (2011) apresentou uma 
metodologia baseada em filtros de Kalman para analise automatica de perturbagoes 
em unidades de geragao de energia com objetivo de fornecer uma estimativa da 
causa das faltas. 

2.3 FERRAMENTAS UTILIZADAS NA EXT R AC AO DE INFORMACOES 

A extragao das informagoes dos registros de oscilografia consiste em urn 
problema tipico de Processamento Digital de Sinais (PDS) e como tal e tratado por 
meio de ferramentas matematicas. Cada ferramenta e mais adequada a urn tipo 
especifico de sinal, portanto o conhecimento dos dados se faz necessario antes da 
escolha por qual ferramenta utilizar, bem como a definigao de quais as informagoes 
a serem extraidas dos dados. A evolugao e adaptagao de metodos matematicos 
aplicados a sinais contmuos no tempo para representar sinais amostrados tern 
ocorrido de maneira rapida e diversas ferramentas surgiram nas ultimas decadas. 
Para citar algumas: filtros digitais, Transformada Discreta Fourier (DFT), 
Transformada Rapida de Fourier (FFT) e a mais recente Transformada Wavelet 
Discreta (DWT) (WEEKS, 2012). 

Diversas publicagoes tern mostrado a aplicagao das tecnicas de PDS nas 
analises de oscilografias. O Quadro 1 relaciona algumas publicagoes de acordo com 
a ferramenta utilizada e o objetivo da aplicagao. 
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Aplicacao 


Ferramenta utilizada 


Referenda 


Identificacao e classificacao de faltas no sistema 
de transmissao da Chesf. 


Transformada de 
wavelet discreta. 


(BRITO, SOUZA e 
COSTA, 201 0) 


Deteccao de disturbios de qualidade de energia. 


Energia dos 
coeficientes wavelets. 


(COSTA et al., 2011) 


Diagnostics em tempo real de disturbios 
transitorios no sistema de transmissao. 


Transformada wavelet. 


(COSTA, 2010) 


Classificacao de disturbios na rede eletrica de 
transmissao de energia. 


Wavelets e redes 
neurais. 


(CANDIDO, 2008) 


Pre-analise e triagem dos registros de 
oscilografias repassados ao ONS pelos agentes 
do sistema de transmissao de energia eletrica. 


Nao informada. 


(GIOVANINI et al., 
2008) 


Analise automatica de disturbios de qualidade de 
energia eletrica. 


Transformada de 
Fourier. 


(KEZUNOVIC e 
LIAO, 2001) 


Monitoracao e avaliacao do estado de 
transformadores de potencial capacitivo. 


Transformada rapida 
de Fourier. 


(LACHMAN, 2007) 


Avaliacao do desempenho de reles digitais de 
protecao de unidades geradoras. 


Nao informada. 


(LENZ, 2009) 


Identificacao e classificacao de faltas em sistemas 
de transmissao de energia eletrica. 


Transformada wavelet 
e rede neural 
probabilfstica. 


(MENESES e 
ROCHA, 2009) 


Analise automatica de perturbacoes em unidades 
geradoras para arquivar automaticamente 
registros que nao necessitem de analise 
detalhada. 


Sistemas especialistas 
e filtro de Kalman. 


(MORETO, 2011) 


Monitoramento e estimativa dos tempos de 
operacao de disjuntores em subestacoes de 
transmissao de energia eletrica. 


Transformada de 
Fourier e transformada 
wavelet. 


(SILVA, 2005) 


Localizacao de faltas em linhas de transmissao. 


Transformada wavelet. 


(SILVA, 

OLESKOVICZ e 
COURY, 2005) 


Deteccao dos instantes de mudangas abruptas 
nos registros de oscilografias da rede de energia 
da Africa do Sul. 


Filtros de Kalman. 


(UKILeZIVANOVIC, 
2007) 



Quadro 1 : Publicacoes: objetivo da aplicacao versus ferramenta utilizada. 
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3. FUNDAMENTOS DA TRANSFORMADA WAVELET 

3.1 POR QUE WAVELETS 

O emprego da matematica para descrever e explicar fenomenos fisicos 
nao e nenhuma novidade. No passado diversos fisicos tiveram que desenvolver 
teorias matematicas para representar fenomenos por eles estudados. No seculo XVII 
o fisico ingles Isaac Newton, sentiu necessidade de ferramentas matematicas 
especificas para exprimir relagoes entre grandezas mecanicas. Ao estudar o que 
chamava de "fluxoes das variaveis" Newton definiu o conceito de derivadas e 
integrals, a base do calculo infinitesimal amplamente utilizado atualmente. 
(CARVALHO e D'OTAVIANO, 2006). 

Nos dias atuais a matematica nao so descreve os fenomenos fisicos, mas 
possibilita a solugao de problemas do mundo real. Atraves do processamento 
computacional solucionam-se problemas que, em alguns casos, necessitariam anos 
para ser solucionados de maneira analitica, alem disso, e possivel obter solugoes 
para problemas que sequer possuem solugao analitica. E quando se trata de 
problemas do mundo real recai-se diretamente sobre o Processamento Digital de 
Sinais (PDS), area da engenharia que tern evoluido com metodos cada vez mais 
complexos e sempre embasados em teorias oriundas da matematica (WEEKS, 
2012). 

No seculo XIX o matematico e fisico frances Jean Baptiste Joseph Fourier 
publicou seu estudo sobre ondas e calor no qual propoem que urn sinal periodico 
qualquer pode ser representado por uma unica expressao analitica. Ele acabava de 
char uma das ferramentas mais utilizadas na analise de sinais e que leva seu nome, 
a "Analise de Fourier" ou no metodo matematico Transformada de Fourier (FT - 
Fourier Transform) (OLIVEIRA, 2007). Aplicando-se a FT a urn sinal optico, por 
exemplo, obtem-se a sua decomposigao nas cores fundamentals do arco iris, onde 
cada cor corresponde a uma frequencia. Se uma das frequencias possuir maior 
amplitude que as demais o feixe de luz tera a cor correspondente a esta frequencia 
predominante. Este e urn exemplo basico da analise de sinais que pode ser 
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desenvolvida a partir do metodo matematico proposto por Fourier e que na 
matematica e representado pela integral da equagao (1). 



fit) : sinal contmuo no tempo t; 
w. frequencia angular em rad/s; 
F(w): sinal contmuo na frequencia w. 
Fonte: (OLIVEIRA, 2007) 

A FT tornou-se uma ferramenta amplamente utilizada em diferentes 



campos como optica, fisica, telecomunicagoes, computagao, eletronica e eletrica. 
Essencialmente ela consiste em uma transformagao linear que troca a base do 
tempo para a base da frequencia de urn sinal. Isso permite olhar, analisar e 
manipular informagoes de frequencia, o que em determinadas aplicagoes e mais 
adequado do que o tempo (MIX e OLEJNICZAK, 2003). 



Ao aplicar a FT em urn sinal contmuo no tempo, e possivel visualizar as 



senoides fundamentals que formam este sinal. A Figura 8 mostra urn exemplo de 
sinal que e composto por varias senoides de diferentes frequencias e amplitudes. 




— oo 



Onde: 



(1) 
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Figura 8: a) Sinai continue* no tempo, b) Decomposicao feita pela FT em sinais de 
frequencias fundamentais. 
Fonte: Autoria propria. 

Um sistema de som permite alterar amplitudes de bandas de frequencias 
especificas em uma musica, realizando sua mixagem, eliminando ruidos 
indesejados ou criando uma nova versao da musica original. Esse tipo de aplicagao 
pode ser feito em um arquivo de audio digital com dados discretos, por exemplo, no 
formato mp3, utilizando a Transformada de Fourier Discreta (DFT - Discrete Fourier 
Transform) que pode ser expressa na forma da equagao (2) (WEEKS, 201 2). 



N-l 



X[m] = V x[n\e-i 2nnrn ' N 



71=0 



Onde: 

x[n]: sinal discreto com (iV - 1) amostras; 

X[m] \ sinal discreto transformado onde m = 0, ...,N - 1; 



(2) 



Fonte: (WEEKS, 2012) 
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Este tipo de analise em sinais discretos e exemplificado na Figura 9 que, 
atraves da aplicagao da DFT, mostra as frequencias que compoem um sinal discreto 
no tempo. 
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Figura 9: a) Sinal discreto no tempo, b) Magnitude dos fasores de frequencias distintas 
que compoem o sinal discreto, obtidas a partir da DFT. 
Fonte: Autoria propria. 

A FT mostra informagoes sobre quais frequencias formam um sinal, 
representando uma media de todo o sinal desde seu inicio ate seu final, ou seja, do 
tempo -°° ao +°° como pode ser visto na integral da equagao (1). A contribuigao de 
cada harmonica do sinal completo fica explicita. Porem a informagao correspondente 
ao instante de tempo em que as frequencias ocorrem e perdida. Isto fica mais 
evidente quando aplicamos a FT em um sinal nao estacionario como o da Figura 10. 
Na Figura 10 b) ve-se apenas a informagao das frequencias que compoem o sinal 
da Figura 1 0 a) mas nao o instante de tempo em que elas ocorrem. 
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Figura 10: a) Sinai de registro de oscilografia. b) Frequencias que compdem o sinal, obtidas a 

partir da FT. 

Fonte: Autoria propria. 

A grande desvantagem da FT para aplicagao em analise de oscilografias 
e o fato de ela nao ser capaz de fornecer informagoes sobre o tempo, apenas 
frequencias. Sendo assim ela e util para analise de sinais estacionarios, ou seja, que 
apresentem urn comportamento repetitivo ao longo tempo, como os sinais da Figura 
8 e da Figura 9 (OLIVEIRA, 2007). 

Para analisar sinais nao estacionarios Denis Gabor, em 1946, propos uma 
adaptagao da FT. Esta adaptagao consiste em dividir, o sinal sob analise, em 
intervalos menores e aplicar a FT a cada intervalo ou "janela" de tempo do sinal 
original, definindo assim a Transformada de Fourier de Tempo Curto (STFT - Short 
Time Fourier Transform) como mostra a Figura 1 1 . 
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Figura 11 : Transformada de Fourier de Tempo Curto (STFT). 
Fonte: Adaptado de (DELMONT, 2003). 
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Para cada janela da Figura 1 1 e aplicada a FT ao sinal e obtem-se como 



resultado a informagao das frequencias que o compoem gerando o grafico 
bidimensional tempo x frequencia. A ST FT consiste em uma transformada local, que 
analisa o sinal no tempo f(t) atraves de uma "janela" W(t), centrada em urn instante 
de tempo t e limitada ha urn tempo inicial e final. A representagao matematica da 
ST FT e dada pela equagao (3) (OLIVEIRA, 2007). 



W(t - t): fungao janela centrada em t e de largura t; 
STFT(w,t): sinal contmuo na frequencia w para cada 
janela. 

Fonte: (OLIVEIRA, 2007) 
A Figura 12 c) mostra o piano tempo versus frequencia, de urn registro de 



oscilografia com perturbagao, obtido pela STFT com uma largura de 0,1333 s para a 
janela de tempo. A Figura 12 a) mostra em detalhe a perturbagao no sinal registrado 
que e mostrado na Figura 12 b). Analisando a Figura 12 c) e possivel identificar que 
entre os instantes 2,1 s e 2,2 s surgem componentes de frequencias maiores que 
nao ocorrem no restante do sinal. Isso pode ser visto na Figura 12 c) atraves da cor, 
cores mais "quentes" representam maior amplitude das componentes de frequencia. 
Atraves do detalhe da Figura 12 a) percebe-se que o intervalo mais exato da 
perturbagao esta entre 2,025 s e 2,045 s. Se a largura t da janela de tempo for 
reduzida podera informar com mais precisao o instante da perturbagao. Porem, ha 
uma dificuldade na escolha da largura da janela, uma vez que, normalmente, nao se 
conhece a-priori o tipo de perturbagao que podera incidir no sinal. 



A STFT relaciona informagoes de frequencias com o tempo, minimizando 



a deficiencia da FT. Verifica-se que a STFT fornece mais informagoes, porem a 
resolugao das informagoes sobre as frequencias e a mesma em todo o intervalo de 
tempo analisado e determinada pela largura da janela escolhida para a STFT. 
(COSTA, 2010). 
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Onde: 



(3) 
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Figura 12: a) Trecho de registro de oscilografia com perturbacao. b) registro de oscilografia de 
curta duracao. c) Piano tempo versus frequencia, obtido pela STFT. 
Fonte: Autoria propria. 



O sinal analisado na Figura 12 e urn sinal tipico dos registros de 
oscilografia de sistemas eletricos. Sinais com caractensticas semelhantes tambem 
ocorrem em registros de sismografos, os quais foram estudados pelo geofisico Jean 
Morlet na decada de oitenta, orientado por Alex Grossman (fisico teorico) na Franga. 
Nestes sinais as componentes de frequencias mudam rapidamente ao longo do 
tempo podendo surgir a qualquer momento. Nestes casos a STFT nao e suficiente 
para extrair informagoes precisas dos sinais. Nos intervalos que contem frequencias 
mais altas uma maior aproximagao e conveniente, necessitando-se variar o tamanho 
da janela para determinar com exatidao o tempo ou a frequencia. Assim Morlet 
propos uma nova transformada, denominada Transformada Wavelet (WT - Wavelet 
Transform), na qual a largura da janela e variavel, permitindo que altas frequencias 
possam ser localizadas com uma maior precisao (DELMONT, 2003). 
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Oliveira (2007) exemplifica com a comparagao: "a analise atraves da WT 
permite visualizar tanto a floresta como as arvores", o que pode ser percebido na 
Figura 13. Na Figura 13 b) ve-se que a janela de tempo que localiza as frequencias 
mais altas e menor do que na Figura 13 a), ou seja, a resolugao para as frequencias 
mais altas e melhor na WT do que na STFT. 
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Figura 13: Resolugao no piano tempo-frequencia a) Obtido pela STFT b) Obtido pela WT. 



Fonte: Adaptado de (OLIVEIRA, 2007). 



3.2 A TRANSFORMADA WAVELET CONTINUA - CWT 



Matematicamente a Transformada Wavelet Contmua (CWT - Continuos 
Wavelet Transform) pode ser expressa pela equagao (4). 



F(a,b) = f(t)V aJ ,(t)dt 

J — 00 

Onde: 

/(t): sinal contmuo no tempo t; 

^a.bi^y wavelet mae, base da transformada. 



(4) 



Fonte: (DELMONT, 2003) 



Em (4) os parametros a e b variam continuamente no conjunto dos reais 
E, sendo a fungao W a b chamada de wavelet mae e definida conforme a equagao (5) 
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1 (t - b\ 

^(t) = -=y — 



Onde: 

a: parametro de escala ou dilatagao; 
b: parametro de tempo ou translagao. 



(5) 



Fonte: (DELMONT, 2003) 



Atraves das equagoes (4) e (5) verifica-se que a CWT depende de dois 
parametros, a e b, que correspondem, respectivamente, as informagoes de escala e 
de tempo, ou dilatagao e translagao. Desta forma urn sinal unidimensional f(t) e 
transformado em urn sinal bidimensional F(a,b) permitindo uma analise tempo 
versus frequencia assim como na STFT, porem com a resolugao aumentada quanto 
maior forem as frequencias, fornecendo uma especie de zoom na escala de tempo 
das mais altas frequencias. (DELMONT, 2003). 

Comparando as equagoes (1) e (4) percebe-se que a CWT representa 
uma generalizagao da FT. Tomando a equagao (4) e utilizando como wavelet mae a 
fungao exponencial complexa obtem-se a equagao (1), da FT, como urn caso 
particular da CWT. A CWT e uma transformada linear e sua base e a wavelet mae 
(OLIVEIRA, 2007). 



3.2.1 Farmliasde wavelets 



Diversas fungoes podem ser escolhidas como base ou wavelet mae para 
a CWT, porem devem ter comportamento oscilatorio e de curta duragao. A origem do 
termo wavelet se deu pela uniao dos termos: "oscilatorio" (onda) e "curta duragao" 
(ondinha). Dai urn dos criterios para verificar se uma fungao pode ser utilizada como 
wavelet e provar que ela e oscilatoria, ou seja, que seu valor medio e nulo. Isto pode 
ser feito atraves da integral da equagao (6). 



+ oo 



<F(t) = 0 



(6) 



Fonte: (OLIVEIRA, 2007) 
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A wavelet mae tambem deve possuir energia limitada, ou seja, deve ser 
verificada a condigao da equagao (7). 



I 



+00 



<F 2 (t) < 00 



(7) 



Fonte: (OLIVEIRA, 2007) 



Diversas fungoes tern sido estudadas e propostas como wavelets. A 
wavelet mais simples e a Wavelet de Haar, que consiste em uma fungao contmua 
por partes e pode ser empregada na analise de sinais com mudangas abruptas. A 
equagao (8) define a Wavelet de Haar e a Figura 14 a ilustra. 



( +1 para (0 < t 0,5) 

W (haar)^ = Ij para ^ < t < X) 

v 0 caso contrkrio 



(8) 



Fonte: (DELMONT, 2003) 




Figura 14: Wavelet de Haar. 
Fonte: Gerada com MATLAB" 



A Wavelet Chapeu Mexicano e definida pela equagao (9) e sua forma, ja 
anunciada pelo nome, pode ser vista na Figura 15. 

~< ^ 2(t 2 - l)e- t2 ' 2 I 

Fonte: (OLIVEIRA, 2007) 
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Figura 15: Wavelet Chapeu Mexicano. 
Fonte: Gerada com MATLAB®. 



A Wavelet de Morlet foi uma das primeiras propostas com a finalidade de 
analise de sinais e se originou dos estudos de Morlet sobre geofisica e analise de 
sinais sismicos. A Wavelet de Morlet esta definida na equagao (10) e sua forma 
exemplificada pela Figura 16. 



if/(morl)^ _ 



7T 



1/4 



e -t 2 /2 e -jw 0 t 



(10) 



Fonte: (OLIVEIRA, 2007) 




Figura 16: Wavelet de Morlet. 
Fonte: Gerada com MATLAB 4 
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A familia de Wavelets de Daubechies foi descrita, no fim da decada de 80, 
pela pesquisadora Ingrid Daubechies e e uma das mais utilizadas na analise de 
oscilografias (DELMONT, 2003). As Wavelets de Daubechies nao sao descritas por 
fungoes explicitas, mas sim criadas por metodos iterativos (OLIVEIRA, 2007). A 
Figura 17 mostra a familia de Wavelets de Daubechies que, em homenagem a sua 
autora, sao nomeadas por "db" seguidas por urn numero "N" que corresponde a 
ordem da wavelet. 
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Figura 17: Wavelets de Daubechies. 
Fonte: Gerada com MATLAB®. 
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A gama de possibilidades para wavelets e infinita e muitas familias ja 
foram propostas para diferentes finalidades. O Quadro 2 cita algumas das familias 
de wavelets mais conhecidas e difundidas na literatura atual. 



Sigla de representacao da familia 


Nome da familia de wavelet 


haar 


Wavelet de Haar 


db 


Wavelet de Daubechies 


sym 


Wavelet Sym lets 


coif 


Wavelet Coif lets 


bior 


Wavelet Biortogonal 


gaus 


Wavelet Gaussiana 


mexh 


Wavelet Mexicana 


meyr 


Wavelet de Meyer 


morl 


Wavelet de Morlet 



Quadro 2: Familias de wavelets. 
Fonte: Adaptado de (MISITI et al., 2014). 



3.3 A TRANSFORMADA l/V/WEZ-ET DISCRETA - DWT 



A WT em sua forma continua (CWT) e util essencialmente para estudos 
teoricos, dedugoes e verificagoes de suas propriedades. Para aplicagoes praticas a 
discretizagao da WT e mais eficiente e permite a sua implementagao computacional. 

Na Transformada Wavelet Discreta (DWT - Discrete Wavelet Transform) 
os parametros de dilatagao e translagao nao variam continuamente, como no caso 
da transformada continua, mas sim discretamente. Em certas aplicagoes, incluindo 
aquelas em analise de sinal, podem-se restringir os valores dos parametros a e b 
(da equagao 5) a uma grade discreta, fixando urn passo de dilatagao a 0 e urn passo 
de translagao b 0 . A equagao 1 1 mostra a forma discreta da transformada em 
comparagao com a forma continua, primeiramente discretizando apenas os 
parametros de escala e translagao (a e b) e depois discretizando tambem a variavel 
independente (tempo). 
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1 (t — b\ 
V aib (t) = — V — — 
Va V a ) 



1 /t - n£ 0 a 0 ? 
<F m , n (t)=— ==W 



V a 0 m J 



1 ft - b\ r , 1 //c - n/j 0 a 0 m \ 
f a ,,(t) = — f — - y fc = y( 2-2- 



Onde: 



men sao inteiros; 

a 0 > 1: parametro de escala ou dilatagao fixo; 
b 0 : parametro de tempo ou translagao fixo; 
t: variavel independente contmua. 
k: variavel independente discreta. 



Fonte: (OLIVEIRA, 2007) 



Uma maneira eficiente de aplicar esta transformada e atraves de filtros, 
tecnica desenvolvida por Mallat em 1986, onde se tern a decomposigao da wavelet 
implementando a analise multiresolugao. A analise multiresolugao e uma tecnica que 
permite analisar sinais em multiplas bandas de frequencias (CASTANON, 2003). 



3.4 ANALISE MULTIRESOLUQAO - AMR 



A Analise Multiresolugao (AMR), tambem conhecida como algoritmo de 
Mallat, e urn metodo proposto para implementar a DWT. O detalhamento da relagao 
entre a DWT e a AMR extrapola o escopo deste estudo e pode ser consultado em 
(WEEKS, 2012) e (OLIVEIRA, 2007). 

A transformada pelo metodo multiresolugao obtem "aproximagoes" e 
"detalhes" de urn determinado sinal atraves da passagem deste por urn filtro passa- 
alta e outro filtro passa-baixa. O filtro passa-alta produz os coeficientes de detalhe e 
o filtro passa-baixa produz coeficientes de aproximagao da decomposigao do sinal. 
Este processo de decomposigao wavelet por meio de filtros e ilustrado na Figura 18 
(PARRAGA, 2002). 
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Figura 18: Decomposicao wavelet com filtros. 
Fonte: Adaptado de (PARRAGA, 2002). 

O coeficiente de aproximagao e uma representagao de baixa frequencia 
do sinal original e o coeficiente de detalhe e a diferenga entre dois coeficientes de 
aproximagao sucessivos. A aproximagao tern a mesma tendencia geral do sinal 
original e o detalhe mostra as componentes de alta frequencia (SILVA, 2008). 

Os filtros utilizados na AMR tern como caractensticas filtrar e reamostrar o 
sinal. A reamostragem, ou downsample, diminui a taxa de amostragem original do 
sinal, eliminando pontos intercalados. Reamostrar por urn fator dois significa que a 
cada dois pontos urn e eliminado. O resultado e urn sinal com metade do numero de 
amostras do original, onde uma amostra "n" do sinal reamostrado corresponde a 
amostra "2n" do sinal original, porem o instante de tempo associado a essa amostra 
e o mesmo (PARRAGA, 2002). 
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Figura 19: Exemplo de downsample. a) Sinal original, b) Reamostrador. c) Sinal reamostrado. 
Fonte: Adaptado de (MIX e OLEJNICZAK, 2003). 

O downsample pode ser realizado multiplicando urn sinal S[n] por urn 
sinal P[n] tal que P[n] seja formado por uma sequencia de zeros e uns, o que resulta 
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em um sinal Ss[n] comprimido e sem perder informagoes como exemplifica a Figura 
19 (MIX e OLEJNICZAK, 2003). 

O processo de decomposigao wavelet da analise multiresolugao 
corresponde a uma decomposigao em sub-bandas, que pode ser implementada com 
um banco de filtros que realizam uma sucessiva decomposigao e reamostragem nos 
coeficientes de aproximagao. Atraves do processo ilustrado na Figura 20, e que 
tambem e chamado de arvore de decomposigao wavelet, obtem-se uma analise 
aproximada de todo o espectro de frequencias (SILVA, 2008). 
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Figura 20: Decomposigao com banco de filtros. 
Fonte: Autoria propria. 

Em termos de algebra linear a decomposigao com banco de filtros cria 
uma sequencia de espagos vetoriais aninhados ou subconjuntos, de tal maneira que 
cada etapa da filtragem corresponde a um zoom no sinal original, como exemplifica 
a Figura 21 . 



V 0 3 V x =J V 2 => ••• ^ V m 



Para m=2: 




Figura 21 : Decomposigao por banco de filtros na forma de espagos vetoriais. 
Fonte: Adaptado de (PARRAGA, 2002). 
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3.5 APLICAQAO 

A fim de exemplificar, uma analise wavelet foi aplicada no mesmo registro 
de oscilografia usado na analise por meio da STFT na Figura 12. Este registro e um 
sinal de tensao registrado durante uma desenergizagao de um gerador, onde a 
tensao decresce lentamente. Na sequencia deste topico sao mostrados alguns 
graficos e figuras que auxiliam na compreensao do metodo. Para realizar a analise 
foi utilizado o Wavelet Toolbox do MATLAB®. 

A wavelet mae empregada para esta analise e a Daubechies de 5- ordem 
(db5) por ser uma das mais empregadas na analise de oscilografias (DELMONT, 
2003). A decomposigao foi aplicada cinco vezes nos coeficientes de aproximagao, 
ou seja, ate o quinto nivel, produzindo cinco coeficientes de detalhes. A Figura 22 
mostra a arvore wavelet da analise realizada. 




Figura 22: Arvore de decomposigao wavelet ate o quinto nivel. 
Fonte: (MISITI et al., 2014). 



Atraves da representagao dos coeficientes de detalhes em escala de 
cores da Figura 23 pode-se visualizar a ocorrencia da escala maxima coincidindo 
com o instante de mudanga abrupta no sinal original. 



44 





' IIIIIIF 

4 - 

3 - 

T 






- 


Q.5 1 1.5 2 


2.5 3 3.5 4 4.5 



Figura 23: Visualizacao da magnitude dos coeficientes de detalhe por escala de cores. 
Fonte: Obtida com Wavelet ToolBox. 



Atraves da visualizagao da magnitude dos coeficientes de detalhes, os 
quais representam altas frequencias, e possivel identificar quais amostras contem 
informagoes de variagoes bruscas no sinal original. A Figura 24 mostra os 
coeficientes de maneira separada. Observando os coeficientes de detalhes percebe- 
se que estes tern predominancia de magnitude quase nula, com excegao para a 
proximidade do instante que ocorre a variagao do sinal original, quando apresentam 
urn pico de magnitude. 

Decomposition at level 5 : s = a5 + d5 * d4 + d3 + d2 * di . 




0 



Figura 24: Visualizagao da magnitude dos coeficientes de detalhe. 
Fonte: Obtida com Wavelet ToolBox. 
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Fica evidente a possibilidade de utilizar as informagoes dos coeficientes 
de detalhes para determinar os instantes de ocorrencias de perturbagoes nos sinais 
de oscilografias, permitindo realizar a segmentagao destes em antes, durante e 
depois de uma perturbagao (UKIL e ZIVANOVIC, 2007). 
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4. ESTUDOS DE CASO 



Para verificar as reais possibilidades de aplicagao da DWT e dos seus 
coeficientes na analise de oscilografias, foram realizados dois estudos de caso com 
registros de ocorrencias reais em unidades geradoras. Os topicos seguintes 
descrevem suscintamente a origem dos registros, mostram uma visualizagao dos 
mesmos e os resultados obtidos com a aplicagao da DWT. 



4.1 CASO 1 : DESLIGAMENTO FORCADO POR CURTO CIRCUITO 



Este caso faz parte da base de oscilografias disponibilizadas para estudo 
por uma concessionaria de geragao de energia do SIN, e foi escolhido por 
representar o desligamento forgado de uma unidade geradora devido a urn curto- 
circuito. 

As caractensticas do sistema sao descritas no Quadro 3 e o ponto de 
ocorrencia da falta foi o para-raios da fase B no lado de alta tensao do transformador 
elevador da unidade. 



Usina termoeletrica movida a carvao 

Unidade geradora 2 

Potencia ativa: 12MW 

Potencia reativa: 5,4 Mvar 

Tensao de operagao: 6 kV 

Frequencia de operagao: 50 Hz 

Gerador conectado em triangulo 

Transformador elevador 6kV para 69kV 

Lado de alta do transformador conectado em estrela 

Quadro 3: Caractensticas da unidade geradora. 

Fonte: Relatorio da ocorrencia disponibilizado para estudos. 

O relatorio disponibilizado pela concessionaria descreve a ocorrencia da 
seguinte forma: 

Desligamento automatico da unidade 2 da UT... provocado pela atuacao da 
protecao diferencial da unidade. Os registros oscilograficos e a sequencia 
de eventos indicam a ocorrencia de urn curto-circuito na fase B do lado de 
alta do transformador da unidade 2, caracterizado por uma sobrecorrente de 
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2696A na fase citada anteriormente. Verificou-se uma subtensao nos 
terminais das unidades associada a sobrecorrentes de cerca de 8000A no 
instante da falta. Os valores pre-falta das correntes das unidades eram de 
aproximadamente 1200A. Apos a eliminacao da falta com a abertura do 
disjuntor da unidade 2, as demais unidades continuaram a operar 
normalmente. 

...houve a explosao do para-raio localizado no lado de alta do transformador 
da unidade 2, o que provocou o curto circuito para terra na fase B. 

A Figura 25 mostra o registro de oscilografia de curta duragao das 
correntes e tensoes da unidade geradora 2, envolvida na falta. Devido a ausencia de 
algumas caractensticas do sistema em que a unidade estava conectada os dados 
sao apresentados em suas unidades primarias e nao estao normalizados em p.u.. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 

Tempo t[s] 



fase A 
fase B 
fase V 



i c 

6.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 

Tempo t[s] 



Figura 25: Visualizacao das correntes e tensoes da unidade geradora 2. 
Fonte: Autoria propria. 




O registro de corrente da fase "B" foi submetido a analise com a DWT, 
com a wavelet mae db5. Os coeficientes de detalhe ate o quinto nfvel de 
decomposigao sao mostrados na Figura 26. 



10 



48 



0 pr 



-10 



CM n 
"O U 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



x 10 



1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 



500 



1000 



1500 



2000 



2500 



3 0'** 



200 



400 



600 



800 



1000 



1200 



LO 



100 



200 



300 

Amostra n 



400 



500 



600 



Figura 26: DWT da corrente de falta. 
Fonte: Autoria propria. 
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Na Figura 26 os coeficientes estao plotados com a mesma escala de 
amplitude para evidenciar as diferentes caractensticas de cada coeficiente. A 
escolha de qual coeficiente utilizar para determinar o instante da falta tambem passa 
pela escolha da wavelet mae mais adequada. De acordo com Costa (2010) a farmlia 
das wavelets Daubechies tern sido amplamente utilizada na analise de oscilografias 
e dentre elas a db4 e a db6 tern se destacado. 

A Figura 27 mostra os coeficientes d2 e d3 em uma escala mais 
apropriada as suas amplitudes. Com urn zoom na escala das amostras evidenciando 
a ocorrencia de urn pico no valor dos coeficientes nas amostras associadas ao 
instante da falta. 
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Figura 27: Coeficientes de detalhe d2 e d3. 
Fonte: Autoria propria. 

Para identificar o instante de tempo da falta pode-se relacionar a amostra 
de maximo valor com o tempo associado a esta amostra, fornecido pela oscilografia. 
Desta forma para o coeficiente d2 o instante de tempo corresponde ao indice t [4n ] e 
para o coeficiente d3 a t[8n]- Porem como os valores dos coeficientes podem assumir 
valores negativos e vantajoso calcular a sua energia. Este calculo e feito atraves da 
integral discreta das amostras dos coeficientes, conforme mostra a equagao 12. 
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Energia do sinal original: 

ife=n 

E 



ZK — ll 
k=l 



Energia dos coeficientes de detalhe: 

k=n/2i , 



ZK = Il 

Onde: 

/: sinal original; 

n: numero de amostras do sinal original; 

k: mdice da amostra atual; 

df. coeficiente de detalhe; 

j[k]\ mdice do coeficiente de detalhe; 



(12) 



Fonte: (COSTA, 2010) 



A Figura 28 mostra a comparagao da energia calculada para o coeficiente 
de detalhe d3 em comparagao com os valores assumidos pelo coeficiente para as 
amostras proximas do momento da falta. 
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Figura 28: Energia do coeficiente de detalhe d3. 
Fonte: Autoria propria. 

A energia dos coeficientes aumenta de forma brusca no momento da 
ocorrencia de transitorios no sinal original e, de acordo com Costa (2010), a 
identificagao deste pico de energia e uma forma apropriada para a determinagao do 
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instante de ocorrencia do transitorio e consequentemente da falta. A partir da 
determinagao do instante que ocorreu a falta pode-se proceder com a segmentagao 
do registro em pre-falta e pos-falta extraindo caractensticas uteis no diagnostico. 

Para detectar o ponto de transitorio, ou momento da falta, atraves da 
energia dos coeficientes de detalhe basta calcular a diferenga entre os valores de 
energia para amostras sucessivas do coeficiente, ou ainda, adotar uma janela de ate 
quatro coeficientes. Assim pode-se estabelecer urn limiar de aumento de energia e 
quando este for ultrapassado determina-se em qual amostra ocorreu o transitorio. 
Costa (2010) utilizou como limiar, para considerar urn aumento brusco de energia, 
que a energia de uma amostra seja no rmnimo tres vezes maior que o valor de 
energia da quarta amostra anterior, como mostra a equagao 13. A Figura 29 mostra 
o aumento da energia do coeficiente d3 no momento da falta. 



Limiar da diferenga de energia entre amostras consecutivas: 

E(n) — Efn_ 4 ) 



E(n-4) 



> 3 



Onde: 



E (n) : energia da amostra atual; 

E (n _ 4) : energia da quarta amostra anterior; 

(n - 4): janela de quatro amostras; 

Fonte: (COSTA, 2010) 
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Figura 29: Aumento brusco da energia do coeficiente de detalhe d3. 
Fonte: Autoria propria. 



4.2 CASO 2: DESLIGAMENTO POR PERDA DE EXCITAQAO 
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Para testar a aplicagao da analise wavelet em registros de desligamentos 
nao ocasionados por curto circuito utilizou-se um registro de desligamento causado 
por perda de excitagao em uma unidade geradora. Este registro tambem foi 
disponibilizado para estudo por uma concessionaria de geragao de energia do SIN. 

A Figura 30 mostra o registro de oscilografia de curta duragao das 
correntes e tensoes da unidade geradora 





1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 
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Figura 30: Visualizacao das correntes e tensoes da unidade geradora. 
Fonte: Autoria propria. 



O registro de tensao da fase "A" foi submetido a analise com a DWT, com 
a wavelet mae db5. Os coeficientes de detalhe ate o quinto nivel de decomposigao 
sao mostrados na Figura 31 . 
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Figura 31 : DWT da tensao fase A. 
Fonte: Autoria propria. 
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Na Figura 31 os coeficientes estao plotados com a mesma escala de 
amplitude para evidenciar as diferentes caractensticas de cada coeficiente. 

A Figura 32 mostra os coeficientes 62 e d3 em uma escala mais 
apropriada as suas amplitudes. Com urn zoom na escala das amostras evidencia-se 
a ocorrencia de urn pico no valor dos coeficientes nas amostras associadas ao 
instante do desligamento. 
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Figura 32: Coeficientes de detalhe d2 e d3. 
Fonte: Autoria propria. 



A Figura 33 mostra a comparagao da energia calculada para o coeficiente 
de detalhe d3 em comparagao com os valores assumidos pelo coeficiente para as 
amostras proximas do momento do desligamento. 
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Figura 33: Energia do coeficiente de detalhe d3. 
Fonte: Autoria propria. 



A Figura 34 mostra o aumento da energia do coeficiente d3 no momento 
do desligamento. 
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Figura 34: Aumento brusco da energia do coeficiente de detalhe d3. 
Fonte: Autoria propria. 
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5. CONSIDERACOES FINAIS 

A realizagao deste estudo permitiu a revisao de conceitos matematicos e 
de processamento digital de sinais, em especial, adquirir novos conhecimentos 
teoricos sobre a transformada wavelet. A transformada wavelet pode, 
simplificadamente, ser definida como uma generalizagao da transformada de 
Fourier. 

Com a aplicagao da transformada wavelet nos estudos de casos verificou- 
se a eficiencia do metodo para a analise de oscilografias de unidades geradoras. A 
farmlia de wavelets de Daubechies e a mais utilizada em analises de sinais nao 
estacionarios. Os coeficientes de detalhes contem informagoes que permitem 
identificar o instante da ocorrencia de transitorios em registros de oscilografias de 
unidades geradoras de energia eletrica. O calculo da variagao de energia dos 
coeficientes e urn metodo eficiente para localizar os transitorios e pode ser 
implementado atraves de algoritmos. 

Nos testes realizados com a wavelet mae db5 o coeficientes de detalhe 
d3 se mostrou adequado para ser utilizado em algoritmos no calcular da variagao de 
energia entre as amostras. Assim e possivel identificar a amostra em que ocorre 
variagao brusca de energia e determinar o instante de tempo que ocorreu a falta, 
dividindo o registro para extrair caractensticas pre e pos falta. 
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